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Abstract 

Nanoelectromechanical systems (NEMS) are nano-to-micrometer scale mechanical resonators 
coupled to electronic devices of similar dimensions. NEMS show promise for fast, ultrasensitive 
force microscopy and for deepening our understanding of how classical dynamics arises by approx- 
imation to quantum dynamics. This article begins with a survey of NEMS and then describes 
certain aspects of their classical dynamics. In particular, we show that for weak coupling the ac- 
tion of the electronic device on the mechanical resonator can be effectively that of a thermal bath, 
this despite the device being a driven, far-from-equilibrium system. 
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I. INTRODUCTION 



A nanoelectromechanical system (NEMS) consists of a nanometer-to-micrometer (mi- 
cron) scale mechanical resonator that is coupled to an electronic device of comparable di- 
_ QflflflflQ. The ^ _ m a y have a Slm p,e „, sueh a S a 
cantilever (suspended beam clamped at one end) or a bridge (suspended beam clamped at 
both ends) and is fashioned out of materials such as silicon using similar lithographic tech- 
niques to those employed for fabricating integrated circuits. Because of their (sub)micron 
size, the mechanical resonators can vibrate at frequencies ranging from a few megahertz 
(MHz) up to around a gigahertz (GHz) [3,0; we are not normally accustomed to the idea 
of mechanical systems vibrating at such high radio-to-microwave frequencies. 

The coupling to the electronic device can be achieved electrostatically by applying a 
voltage to a metal film deposited on the surface of the mechanical resonator. One example of 
a coupled electronic device is a single electron transistor (SET) shown in figure H3 Electrons 
quantum tunnel one at a time across the transistor from drain electrode to source electrode, 
driven by a drain-source voltage Vd s (we adopt the usual convention where positive charges 
flow from source electrode to drain electrode and hence negative charges flow the opposite 
way). The magnitude of the resulting current depends on another voltage applied to a third, 
gate electrode, called the gate voltage V g . With the metallized mechanical resonator forming 
part of the gate electrode, motion of the former will modulate the gate voltage and hence 
the drain-source tunneling current, which is subsequently amplified and detected. 

With the high frequencies and small inertial masses of the nanomechanical resonators, 
together with the ultrasensitive mechanical displacement detection capabilities of the cou- 
pled electronic devices, NEMS show great promise for metrology. One possible area of 
application is force microscopy, where a cantilever tip is scanned over a surface and the 
cantilever displacements measured as the tip interacts with the surface used to build up a 
force topography map. Of particular interest is the magnetic resonance force microscope 
(MRFM) which employs a ferromagnetic cantilever tip, enabling the mapping of unpaired 
electron and nuclear spin densities at and below the surface |9j]. Recently, single electron 
spin detection sensitivities were achieved 10,[ll|; the potential applications of being able 



to determine chemical identity at the single molecule or atom level are numerous. And 
with the use of smaller, suitably-engineered NEMS MRFM devices, the higher mechanical 
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frequencies might result in faster read-out times at equivalent or better sensitivities. 

Another application is mass-sensing, where the mass of a small particle attaching itself to 
a nanomechanical resonator can be determined from the resulting vibrational frequency-shift 
of the resonator. Recently, attogram ('atto'= 10~ 18 ) detection sensitivities were achieved 3, 
Q|. With the use of suitably-engineered higher frequency NEMS, the detection of individual 
molecules may be possible at single-Dalton sensitivities (1 Dalton = 1.66 x 10~ 27 kg, 1/12 
the mass of a C 12 atom) [l^ |. 

NEMS can be interesting in their own right as nontrivial dynamical systems. Because 
of the small inertial mass of the nanomechanical resonator and the strong electrostatic cou- 
pling to the closely integrated electronic device that can be achieved, individual electrons 
travelling through the latter can give significant displacement 'kicks' to the mechanical res- 
onator. In turn, the motion of the resonator will influence the electron current and so 
on. At cryogenic temperatures, certain electronic devices can behave in a quantum coher- 
ent fashion, existing in a quantum superposition of different charge states as electrons are 
transmitted through the device. Interacting with such a device, the mechanical resonator's 
centre-of-mass may be driven into a quantum state jfj], such as a superposition of separated 
position states. The quantum nature of the coupled electromechanical system will manifest 
itself in certain signatures of the measured current. Nanomechanical resonators comprise 
up to about ten billion atoms, so that by most standards such quantum effects would be 
deemed macroscopic. It is important to appreciate that we are here referring to quantum 
effects in 'dirty real' devices that possess both many electronic and mechanical degrees of 
freedom, and which interact strongly with the surrounding environment consisting of pho- 
tons, phonons, fluctuating (charge) defects in both the mechanical resonator and electronic 
device etc. The experimental and theoretical investigation of such systems will lead to a 
deeper understanding of how classical dynamics emerges as an approximation to quantum 
dynamics; NEMS straddle the microscopic quantum and macroscopic classical worlds. 

In the fir st g eneration of experiments to probe the dynamics of NEMS (see, e.g., Refs. jisl 
Q [ill Q, Q, H- the mechanical resonator component has been found to behave 

classically, as might be expected; the experiments are not yet quite refined enough to observe 
quantum interference effects that are washed-out by the resonator's environment. Despite 
this, the (semi) classical dynamics of NEMS has been found to be nontrivial and worthy 
of investigation. One line of investigation is to identify common features in the classical 
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dynamics of the various NEM devices, so as to bring some degree of coherence to the 
field. Remarkably, under certain conditions of weak coupling and also wide separation of 
mechanical and electronic dynamical time-scales, the electronic device behaves effectively as 
a thermal bath: the mechanical resonator undergoes thermal brownian motion, characterized 
by a damping constant and effective temperature which are determined by the electronic 



parameters of the device |22j, 23]. The fact that the electronic device can be effectively 
replaced by a thermal bath is at first sight rather surprising, given that the voltage-driven 
electron current flowing through the device is a far-from-equilibrium, many electron state. 
The use of well-understood equilibrium concepts in the construction of theoretical models 
of less-well-understood non-equilibrium systems goes back to the early days of statistical 
mechanics and continues to find broad application. For a recent review, see for example 
Ref. Q- 

The outline of this article is as follows: section |H] gives examples of various representative 
NEM devices that are being investigated. Section ITTT1 analyses the classical dynamics of the 
SET-mechanical resonator system, focusing on the effective equilibrium description in the 
weak-coupling regime. Section IIVI describes the effective equilibrium dynamics of some of 
the other NEMS introduced in section |H] Section IS gives concluding remarks. 



II. SURVEY OF NEMS 



In this section we describe several representative NEMS. The explanations as to how they 
work will be qualitative and brief in nature, centering on schematic diagrams or electron 
microscope images of each device. 

Figure ^ shows a schematic diagram of a mechanically compliant tunneling electrode 
taken from Ref. |2£|. The tunneling electrode consists of a cantilever with metal tip which 
is placed close to the surface of a bulk metal counterelectrode. As is common in theoretical 
analyses of NEMS, the cantilever is simply modeled as a harmonic oscillator with some 
effective spring constant and mass. Electrons tunneling across the gap between the metal 
tip and surface cause the cantilever to recoil, while in turn the cantilever's motion affects 
the electron tunneling probability and hence the measured tunneling current. Thus, we 
have a coupled electromechanical system which should display some interesting dynamical 
properties, depending on the applied voltage across the gap, gap size, cantilever's mass, 
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spring constant etc. There have been several 



tunneling electrode and related schemes 
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theoretical investigations of the mechanical 
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32j. In a sense, the 



scanning tunneling microscope [33] is an experimental realisation, since a tunneling electrode 
cannot be made absolutely rigid. However, a device with a deliberately compliant, low-mass 
tunneling electrode, such that the tunneling electrons themselves cause significant recoil of 
the electrode, has yet to be demonstrated. 




FIG. 1: Scheme for a mechanically compliant tunneling electrode j25l |. The cantilever electrode is 
modeled as a harmonic oscillator with some spring constant k c and mass m c . A voltage is applied 
across the gap 2a and the resulting current measured. 



Figure 121 shows a scanning electron microscope (SEM) micrograph of a quantum point 
contact (QPC) displacement detector developed by Andrew Cleland's group at UC Santa 
Barbara j3J]. The detector comprises a suspended beam etched from a single- crystal gallium 
arsenide (GaAs) heterostructure with a fundamental resonance frequency for out-of-plane 
flexural motion of about 1.5 MHz. Within the beam is a thin layer of free electrons, called 

n 

a two-dimensional electron gas, (see, e.g., Ref |3_5j for an elementary review of 2DEG and 
other low-dimensional semiconductor systems) which forms a current when a drain-source 
voltage is applied across the ends of the beam. Applying also a sufficiently negative voltage 
to the metal gate electrodes on the surface of the beam expels the electrons from directly 
beneath the elecrodes so that they can only flow through the narrow, electrostatically defined 
constriction between the point contacts. The constriction can be made narrow enough such 
that it is comparable to the electrons' Fermi wavelength; hence the name 'quantum point 
contact'. Because GaAs is a piezoelectric material, mechanical strain in the flexing beam will 
induce a polarization electric field in the beam. This induced polarization field has the same 
effect as applying a gate voltage, hence modulating the current flowing through the QPC. In 



5 



turn, the fluctuating electric field due to the flowing current will induce a mechanical strain. 
Thus, the QPC displacement detector is another example of a coupled electromechanical 
system. 




FIG. 2: SEM micrograph of QPC displacement detector, showing the suspended beam with elec- 
trodes on the surface 34]. The gate electrodes are labeled T' and '3'. The point contacts forming 
the narrow constriction are located just above label T'. The drain and source electrodes are labeled 
'5' and '2'. The indicated magnetic field is used to actuate the beam; the Lorentz force on the 
electron current causes flexing displacements in the indicated z direction. 



Figures El IH and E] show various realisations of the single electron transistor (SET) dis- 



placement detector. The device in figure El was developed by Cleland's group [HUlo], while 
the device in figure|3]was developed by Keith Schwab's group at the Laboratory for Physical 
Sciences, U Maryland llfil . Figure El shows a scheme being developed by Herre van der 
Zant's group at Delft [36|. The latter device differs from the first two in that the SET island 
is mechanically-compliant instead of the gate electrode. A closely-related transistor device 
integrating a suspended carbon nanotube was recently demonstrated by Paul McEuen's 
group at Cornell jscj. The basic operating principle of the SET displacement detector is 
illustrated in the cartoon (figure EK)- The source and drain aluminum electrodes are electri- 
cally insulated from the island electrode by an oxide layer. However, because the oxide layer 
is very thin, electrons can quantum tunnel across the oxide barriers, from drain-to-island- 
to-source, constituting a current which is measured after subsequent amplification stages. 
In order for an electron to be able to tunnel from, e.g., the drain to the island electrode, 
the total work done by the drain-source and gate voltages must exceed the accompanying 
change in the stored electric field energy due to the redistribution of the charges on the 
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various electrodes (called the 'charging energy'); the difference between the work done and 
the charging energy gives the net energy gained by the tunneling electron, which must be 
positive. Because of the small sizes of the electrodes, their capacitances are small and so 
the charging energy is large. Hence, a large enough drain-source voltage must be applied in 
order to meet this charging energy cost so that electrons can tunnel through the device. If, 
however, the voltage is not large enough to overcome the charging energy required to put 
more than one electron simultaneously on the island, then only one electron can tunnel on 
and off the island at a time (hence the name 'single electron transistor'). 

The charging energy can be offset by the applied gate voltage, so that varying the gate 
voltage modulates the tunneling current. Alternatively, by applying a fixed voltage to the 
mechanically-compliant beam gate electrode, motion of the latter will also modulate the 
tunneling current. Figure 0Jd shows an example from Schwab's group of the SET-detected 
signal of the mechanical beam undergoing thermal brownian motion. The area under the 
peak after subtracting off the white-noise background (which mostly originates from the 
subsequent amplification stages) gives a measured mean-squared displacement (x 2 ) of about 
2 x 1CT 13 m, with a position detection sensitivity set by the white-noise floor of about 
1CT 13 m. To put these numbers into perspective, this SET is able to detect displacements as 
small as one-thousandth the diameter of a hydrogen atom. Furthermore, such sensitivities 
are within an order of magnitude of the quantum zero-point displacement uncertainty of the 
mechanical beam. 

Using the classical equipartition of mechanical energy for a harmonic oscillator, (E) = 
muj 2 {x 2 ) = k-^T, it is possible to directly determine the temperature of the mechanical beam 
from the measured mean-squared displacement, resonant frequency and effective motional 
mass (which can be estimated from the beam dimensions and mass density). The data in 
figure 0Jd corresponds to a beam temperature of about 60 mK. This value sets a record for 
the lowest directly measured temperature of a nanomechanical resonator. The actual, base 
temperature of the refrigerator is about 35 mK, suggesting that there is some local heating 
of the beam. Most likely, the energy released by the tunneling electrons in the SET is not 
dissipating quickly enough. 

In addition to heating the mechanical beam resonator, the tunneling electrons will di- 
rectly influence the motion of the resonator: the fluctuating island voltage due to electrons 
tunneling on and off the island will produce a back reaction force on the electrostatically 
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coupled resonator. The magnitude of this force noise increases with applied gate voltage. 
The displacement detection sensitivity also increases, so that there is an optimum gate volt- 
age bias point which is neither too large nor too small. On the other hand, if our interest 
is not in displacement detection, but rather the investigation of the coupled dynamics of 
NEMS, then the gate voltage should be increased beyond this optimum bias point so that 
there is a large back reaction on the resonator . Motion of the resonator then modulates 
the tunneling current and in turn the fluctuating current drives the mechanical resonator. 
Experiments are underway in Schwab's group to explore the coupled dynamics in this back 
reaction-dominated regime. 




beam 



FIG. 3: (a) Cartoon illustrating the operation of the SET displacement detector. The indicated 
charging energy levels represent the energy cost arising from the change in stored electric field 
energy as one or more electrons tunnel onto the island; the discreteness of the levels are not a 
quantum effect, but rather the incremental energy cost to put increasing numbers of electrons on 
the island at the same time, (b) False-coloured SEM micrograph showing the suspended, doubly- 
clamped beam and SET The substrate and beam are fashioned from GaAs (blue regions), and 
the SET and beam gate electrodes are thin layers of aluminum (yellow regions), with aluminum 
oxide forming the tunnel barriers. The beam is located 0.25 microns away from the island electrode. 
The measured fundamental flexural frequency for in-plane motion is about 116 MHz. 

Figure El shows a nanomechanical charge shuttle developed by Dominik Scheible at Ludwig 
Maximilians University, Munich, and Robert Blick at the University of Wisconsin . The 
mechanical shuttle element is a nanopillar fashioned out of silicon with a conducting island 
at the top made out of gold. An ac voltage is applied to the source electrode at a frequency 
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FIG. 4: (a) SEM micrograph showing the mechanical beam and SET wire drain-source wire leads 
and island ^J]. Electrons tunnel one at a time across the electrically- insulating junctions located 
at the corners (J). The beam and surrounding substrate are fashioned from a silicon- nitride (SiN) 
membrane, a commonly-used material for NEMS because of its high strength-to-mass density ratio. 
The beam is coated with a layer of gold, forming the gate electrode, while the SET island and leads 
are made from aluminum. The mechanical resonator is located 0.6 microns away from the island. 
The large, fixed gate electrode at the upper left is used for actuation, (b) SET-detected noise power 
spectrum converted to displacements-squared. The Lorentzian peak due to the mechanical beam's 
thermal brownian motion is clearly seen above the amplifier noise background. The measured 
fundamental frequency for in-plane motion is about 19.7 MHz. 




FIG. 5: Scheme for a SET with mechanically-compliant island formed out of a suspended carbon 
nanotube [26j. As a result of the electrostatic coupling between the voltage-biased gate and island, 
the position of the latter fluctuates as electrons tunnel on and off. 
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close to the fundamental flexural frequency of the pillar. When there is an excess charge on 
the island, the ac voltage exerts a force on the pillar, driving it into mechanical oscillation 
at the ac frequency. If the ac voltage amplitude is sufficiently large, then the shuttle will 
deflect close enough to the source and drain electrodes such that electrons can tunnel between 
the island and the electrodes. The pillar then shuttles charge between the two electrodes, 
producing a current that is detected at the drain electrode. The number of electrons that 
tunnel and the tunnel direction depend on the magnitude and sign of the drive voltage 
at the instant of closest approach to a given electrode, i.e., on the phase lag between the 
oscillating mechanical motion and drive voltage. The phase lag in turn depends on the drive 
frequency relative to the fundamental frequency of the shuttle. Thus, the magnitude and 
direction of the shuttle current can be controlled by varying the drive frequency. The charge 
shuttle bears some resemblance to the SET with mechanically compliant island shown in 
figure However, as the name implies, in contrast to the above-described NEM devices 
mechanical motion of the shuttle is essential for an electron current to flow; the electronic 
and m echanical degrees of & eedo m are inextricably linked in tneir d_ Ret. g ive S 
a review of charge shuttle systems. 




FIG. 6: SEM micrograph of the silicon nanopillar with gold island (I) and source (S) and drain (D) 
electrodes 3|. The gate (G) electrode was not used in the experiment. The current is detected 
from the drain electrode after amplification. The nanopillar has a fundamental flexural frequency 
of 367 MHz. 



In this section, we have introduced several examples of NEM devices. The following 
sections will explore certain aspects of their coupled dynamics, with a central goal being 
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to identify commonalities in their dynamical behaviour. In this respect, it will be useful to 
view a NEM device according to the generic scheme of figure The mechanical resonator, 
which will be simply modeled in its lowest fundamental vibrational mode as a harmonic 
oscillator, is an open system that is coupled to two energy reservoirs. The electronic device 
through which an electron current flows constitutes one of the reservoirs where the current 
exchanges energy with the oscillator via the electromagnetic interaction. All degrees of 
freedom apart from those of the electronic device and the oscillator constitute the second, 
'external' reservoir. These degrees of freedom consist, for example, of higher vibrational 
modes, fluctuating defects etc. within the mechanical resonator and air molecules, photons 
etc. impinging the surface of the resonator. These external degrees of freedom are simply 
modeled as one infinitely large thermal equilibrium bath at some temperature T ext . If the 
oscillator is initially in an excited state, then in the absence of the electronic device it 
will lose energy to this bath at a rate which can be expressed as u>/Q ext , where u> is the 
oscillator frequency and Qext is the quality factor. Nanomechanical resonators are typically 
found to have quality factors between 10 3 — 10 5 jjj, so that a resonator oscillates at least a 
few thousand cycles before its amplitude has substantially decayed from the initial excited 
amplitude. 

We shall pay particular attention to the dynamics of the oscillator system as a result 
of its interactions with these two reservoirs. In this respect, we shall envisage a separate, 
idealized direct probe of the oscillator dynamics allowing perfect measurement of its posi- 
tion and velocity coordinates without affecting its dynamics. In actual NEMS experiments, 
information about the mechanical resonator dynamics is obtained by measuring the elec- 
tron current: the electronic device is after all usually intended as a displacement detector. 
For some theoretical investigations about how the dynamics of the mechanical resonator 
manifests itself in the current signal output, see, e.g., Refs. js^IsS, 39 



III. THE SET-NANOMECHANICAL RESONATOR SYSTEM 

In this section we describe the coupled classical dynamics of the SET-mechanical resonator 
system. The coupled dynamics of other NEMS will be considered in section ITVl Referring to 
the generic scheme of figure we shall in the first instance omit the coupling to the external 
thermal bath (i.e., Q ext —* oo); we are interested in elucidating the 'pure', coupled dynamics 
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FIG. 7: Scheme of a generic NEMS. 

of the oscillator interacting with the SET only. For example, is this system stable in the 
sense that the tunneling electron current removes energy from an initially excited oscillator 
causing damping, or do the tunneling electrons dump energy into the oscillator driving it to 
progressively larger amplitudes? The answer is not obvious a priori. 

A. The master equation 

We seek some form of manageable equation of motion which describes the SET-mechanical 
resonator system. The logical starting point is the time-dependent Schrodinger equation 
with Hamiltonian involving the relevant microscopic degrees of freedom, in particular the 
conduction electrons in the various SET electrodes, and the fundamental oscillator mode 
(recall we are neglecting for the time being the external environment of the oscillator mode). 
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One then proceeds through various stages of approximation to arrive at simpler-to-handle 
classical statistical equations of motion along with some conditions for their validity. How- 
ever, for NEMS at least, this procedure still needs to be properly developed and is essential 
to the future goal to achieve a deeper understanding of how classical dynamics emerges from 
quantum dynamics by approximation in these systems. 

In the absence of such a properly-defined procedure, we shall instead use physical intuition 
to guide us directly towards writing down the classical statistical equations of motion. There 
is a long tradition employing such an approach going back to Boltzmann and his famous 
equation. One advantage is that such equations often accurately capture the statistical 
dynamics as verified by experiment. Another particular advantage in our case is that the 
dynamics in the regime of strong electromechanical coupling is easier to analyze than in 
the full quantum theory. The general disadvantage is that the precise conditions for the 
validity of the classical equations of motion are not available, with the consequence that 
important terms can occasionally be overlooked. This is believed not to be a problem for 
the SET-oscillator classical equations which we shall shortly write down. 

What minimal set of coordinates is required to describe the mechanical resonator inter- 
acting with the SET? Modeled as a single oscillator mode, the resonator's state is described 
by a position x and velocity v coordinate. With respect to the oscillator, the relevant SET 
coordinate is the total number N of excess electrons on its island; depending on N, the 
electrostatic force will pull the oscillator towards the island by different amounts. (If we 
were to also describe the SET drain-source current, we would need an additional counter 
coordinate which keeps track of the number of electrons entering (or leaving) the SET). 
Because tunneling is a random process, the coordinates N, x and v will fluctuate. Thus, the 
equations of motion can either take the form of stochastic differential equations involving the 

n 

random, time- varying functions N(t), x(t) and v (t) (see e.g., Ref. [41| for an introduction on 
how to describe stochastic processes), or we can write down a deterministic equation involv- 
ing the probability density function P/v(x, t>, t). In the latter formulation, PN(x,v,t)dxdv 
is interpreted as the probability at time t of picking out from a large number of identical 
SET-oscillator systems (an ensemble), one system with island number N and with position 
and velocity coordinates in the respective intervals [x, x + dx] and [v, v + dv]. We shall choose 
to express the dynamics in terms of the probability density function, rather than in terms 
of stochastic coordinates. Referring to the circuit diagram for the coupled SET-resonator 
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system (figure [HJ), we have 
dP N 



dt 
dt 



{h n , p N } - (r- L + r+) p N + (r+ + v R ) p N+l 
{H N+1 , p n+1 } - (r+ + r R ) p N+1 + (r- L + r+) p n , 



(i) 

(2) 



where Hn(n + i) is the oscillator Hamiltonian for the resonator in the background electrostatic 
potential of the SET with N(N + 1) electrons on the island, {-, ■} is the Poisson bracket and 
^%R) are ^ e e l ec t r on tunneling rates to the left (+) or to the right (— ) across the left (L) 
or right (R) tunnel junctions. Here, we are assuming that the voltages are chosen such that 
the number of excess electrons on the island fluctuates between N and N + 1 only. 






FIG. 8: Circuit diagram of the SET-mechanical resonator system. The capacitances of the left (L) 
and right (R) tunnel junctions are assumed identical, denoted as Cj. The gate capacitance formed 
by the adjacent resonator and island electrodes is denoted C g . As the resonator flexes, the distance 
between the two electrodes varies, hence changing C g . With the resonator electrode kept fixed at 
the gate voltage V g , the changing C g gives rise to a varying island potential, hence modulating 
the tunneling current. For the indicated drain-source voltage (V^ s ) polarity, the usual direction for 
electron flow is from right to left, governed by the tunnel rates F R and T^. 



Equations 



and (J2J) describing the evolution of the probability distribution 



N(N+1) 



[x,v,t) are commonly called 'master equations'. They are coupled, first-order par- 



tial differential equations in the variables x, v and t. Once an initial probability distribution 
has been specified, then these equations uniquely determine the subsequently evolving prob- 
ability distribution. For example, the distribution at some initial time to might take the form 
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Pn{x, f, to) = S(x)S(v) and Pn+i(x, v , to) = 0, meaning that at time t all SET-oscillator 
systems within the ensemble are in the same state, with the oscillator's position and velocity 
being x = and v = 0, respectively, and the SET island number being N. If the tunneling 
rate terms were absent in equations ([TJ and (J2J) then the evolving probability distribution 
would remain a delta function, describing deterministic harmonic oscillator motion with the 
island number remaining fixed at N. Choosing the origin of the coordinate x to coincide 
with the equilibrium position of the oscillator when there are N electrons on the island, the 
oscillator Hamiltonians take the form 



Hn = — — I — muj 2 x 2 (3) 
2m 2 w 

p 2 1 

fljv+i = - h -mu 2 (x -x ) 2 , (4) 

2m 2 

where Xq is the distance between equilibrium positions of the oscillator with iV and iV + 
1 electrons on the island. With the tunneling rate terms present, however, the evolving 
probability distribution begins to spread. The fluctuating electron island number due to 
electrons tunneling on and off the island causes the deterministic evolution of the oscillator 
to be interrupted at random times by a sudden shift ±xo in the origin of the harmonic 
potential experienced by the oscillator. 

The placings and the signs in front of the various tunneling rate terms in the master 
equations (|H) and (j2J are easily understood from figure |H1 For example, in equation (JTJ), 
the first two rates and correspond to tunneling onto the island and so decrease the 
likelihood that the island number remains at N; hence the 'minus' sign. On the other hand, 
the second two rates and correspond to tunneling off the island and so increase the 
likelihood that the island number becomes N; hence the 'plus' sign. The tunneling rates 
can be derived using Fermi's Golden Rule j3] and take the form 

r ± ... 1 ^L(R) 

^(K) e*Rj 1 _ e ~E ± L(R) /^ W 

where Rj is the effective tunnel junction resistance (assumed the same for each junction), 
T e is the temperature of the source, drain and island electron reservoirs (assumed the same) 
and E L , R ^ is the energy gained by a single electron tunneling to the left (+) or to the right 
(— ) across the left (L) or right (R) junction. When the electron temperature T e is small 
compared to the single electron charging energy, i.e, k^T e <C e 2 /2C*s (where total SET 
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capacitance is Cs = 2Cj + C g ), then equation (j3J) becomes approximately 

T t(R) = -J^ E L{R) e ( E L{R)^ ( 6 ) 

where G(-) is the Heaviside step function. Thus, a given tunnel rate is only non-negligible 
provided the associated E is positive. For micron-scale SETs, ~ 1 fF (= 10~ 15 F), giving 
T e C 1 K. In most experiments, this condition is satisfied and so from now on we will use 
expression © for the rates. For the indicated drain-source voltage polarity in figure |H1 rates 
rjj and describing tunneling to the left are non- negligible, while the rates and 
for tunneling to the right are exponentially suppressed. Note, however, that through the 
£"s, the tunnel rates also depend on the oscillator's position x; it is possible that, for large 
enough amplitude displacements and applied gate voltages, a given E changes sign hence 
switching on or off the corresponding tunnel process and changing the tunneling direction. 

The above master equation pair (JJJ and (J2J) is clearly classical. The probability distri- 
bution characterizes the classical statistical uncertainty in ones knowledge of the precise 
coordinates x, v of the oscillator and island electron number. The equations do not en- 
tertain the possibility of quantum superpositions between different island number states, 
different oscillator position/velocity states, or entangled oscillator-island states. Quantum 
mechanics only enters in the determination of the effective tunnel junction resistances Rj 
appearing in the rate expressions (JHJ), essentially providing the randomness of the transi- 
tion rates. Refering to the discussion at the beginning of this section, the classical master 
equation can in principle be derived from the Schrodinger equation describing the time- 
evolution of the density matrix characterising the quantum state of the harmonic oscillator 
and conduction electrons in the SET electrodes. One step in the derivation is to 'trace-out' 
the many electron state, leaving in the SET state description just the total island number. 
This step introduces irreversibility into the equations, as manifested by the tunneling rate 
terms. The two key approximation steps out of which a classical description emerges are 
(1) assume a large enough SET island such that electron energy level quantization can be 
neglected; (2) assume large enough tunnel junction resistances such that the electron tun- 
nels incoherently between the drain/source and island electrodes. Another approximation 
step that is made is to assume that the timescale for polarization charges on the electrode 
surfaces to re-equilibrate in response to a tunneling event is negligible as compared with the 
characteristic time between tunnel events given by r tunne i = eRj /V& a . This is reflected in the 
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fact that the rate of change of the probability densities at time t in the master equation are 
governed by transition rates that are weighted by probability densities at the same time t 
and not at earlier times. This is called the 'Markovian' approximation. With the oscillator 
excluded, the above master equation arising from the various just- described approximations 
often goes by the name of the SET 'orthodox model' (see, e.g., chapter 4 of Ref. j3| for an 
accessible discussion) . 



B. The steady state solution 

In solving the master equations (JIJ and (J2J), we will first determine the steady state be- 
haviour for t — > oo. Assuming sufficiently large Vd s and sufficiently weak coupling between 
the oscillator and SET such that the energy terms and E£ are always positive, then 
we need only consider the rates and as defined in equations (jBJ) with the step func- 
tion omitted. It is convenient to express the master equations in terms of dimensionless 
coordinates, since in dimensionless form the essential parameters governing the dynamics 
are more clearly expressed. Expressing the time coordinate in tunneling time units r tunnol , 
the position coordinate in shift units xo and the velocity coordinate in units xo/r tunne \, the 
master equations take the form 

P N = e 2 *^ -v^ + EtP N+1 - E%P N (7) 
P N+1 = e\x - 1)^1 _ v ^± _ EiP N+l + E+P N , (8) 

where the dimensionless parameter e = u;r tunne i characterizes the separation between the 
oscillator and SET dynamics timescales. The dimensionless energy terms are obtained by 
dividing E± {R) by eV ds : 

E+ = --L-(N g -N- 1/2) - kN + 1/2 - kx (9) 

Er = +tAt( N 9 ~ N ~ 1/2) + kN + 1/2 + kx, (10) 

where N g = C g V g /e is the polarization charge induced by the gate voltage and the di- 
mensionless parameter k = muj 2 x^/ {eV^) characterizes the coupling strength between the 
oscillator and the SET. The shift coordinate has the explicit form x$ = —eN g / (CsmLJ 2 d) , 
where d is resonator-island electrode gap, so that the coupling strength can be controlled by 
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varying the gate voltage V g ; in particular, k depends quadratically on V g . As we shall see, 
parameters e and k are key to describing the coupled SET-oscillator dynamics. 

From (0) and (|BJ), we can derive equations for the various moments (x n v m ) n(n+i) of the 
probability distribution, where 



(x n v m ) N{N+1) = j dx j dv x n v m P N(N+1) (x,v,t). 

Solving for the moments is a more manageable task than trying to solve for the full proba- 
bility distribution all at once. The equations for the moments are 

(x n v m ) N = -me 2 (x n+1 v m - 1 ) N + n(x n - 1 v m+l ) N + E L (x n v m ) N+l 

-E R (x n v m ) N - k ((x n+l v m ) N + (x n+1 v m ) N+1 ) (11) 

(x"v™) N+1 = -me 2 ((i" +1 . m -% + i - (xV- + «(i n -^ m+ % +1 
-E L (x n v m ) N+1 + E R (x n v m ) N 

+k ((x n+1 v m ) N + (x n+1 v m ) N+1 ) , (12) 

where the K-dependent coupling terms have been pulled outside the energy E terms and 
the '~' and '+' symbols on the latter have been dropped for notational convenience. If 
the SET damps the oscillator, then in the limit t — > oo the various moments approach 
constant values. Thus, we seek a possible solution to (fTTj) and (fT2j) with {x n v m ) n(n+i) = 0. 
Fortunately, these equations can be solved exactly for the time-independent moments. For 
n + m = 0, 1, we find 

(x) N+ i = (P) N+ i = 1 - (P) N = (13) 

1 — K 

and (x)n = (v)n = (v)n+i = 0. For n + m = 2, we find for the variances 2s| 

6x 2 = (x 2 ) - (x) 2 = ^(P) N (P) N+1 (14) 

Sv 2 = (v 2 ) = (1 - K )^(P) N (P) N+1 (15) 

m 

and (xv) at = (xv) n+i = 0, where the averaged probabilities are defined as 

(P)n(n+i) = J du J dv P N{N+1) (x,v,t) 

and we have returned to the original, dimensionful coordinates. 

What can we learn from these solutions to the first few moments? For k < 1, the very 
existence of the solutions is evidence that the SET damps the oscillator, bringing it to a 
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steady state. For k > 1 on the other hand, the solutions are ill-defined (note that 5v 2 , 
Pn+i < if k > 1) suggesting that the coupled SET-oscillator system may be unstable in 
this regime. It is not possible to say anything more definite than this about the behaviour 
in the large k regime, however, since the approximations to the master equation made above 
(in particular the omission of the step functions and the tunneling rates to the right) break 
down when k is not small. Note from (fTlj) and (|15|) that the position and velocity variances 
are related as 

uj 2 r 8x 2 = 5v 2 = ^(P) n (P) n+1 , (16) 
m 

where the renormalized oscillator frequency is ujr = y/1 — k uj. Now recall that for a 
classical damped harmonic oscillator in contact with a thermal bath at temperature T, we 
have equipartition of energy: ^m5v 2 = ^mco 2 5x 2 = \k-oT. Comparing with equation (JH3), 
this suggests that we can assign an effective temperature T SET to the SET given by 

A^bTset = eV ds (P) N (P) N+1 . (17) 

Further evidence that the SET behaves effectively as a thermal bath comes from examining 

the higher moments with n + m — 3, 4 For k <C 1, it is found that the higher moments 

can be approximately decomposed into products of the lowest moments with n + m = 1, 2, 
i.e., (x), (x 2 ), and (v 2 ). This means that the steady state probability distribution P(x,v) [= 
Pn(x, v) + P/v + i(x, v)} for the oscillator is given by a Gaussian to a good approximation, as 
is the Maxwell-Boltzmann distribution describing a classical harmonic oscillator in contact 
with a thermal bath: 

P ( XiV ) = ^^exp (--^- [uj\x - {x)f + v 2 }). (18) 

Figure 01 shows an example steady-state probability distribution for P(x) [= J dv P(x, v)] 
obtained by numerically solving the full master equations (0) and with parameter choices 
e = 0.3 and k = 0.1. Also plotted is the Gaussian distribution obtained after integrating 
(fTS|) over the velocity coordinate, with the steady state average position coordinate given 
by (x) = xo(Pn+i) [the dimensionful form of equation (JT3J)] and the temperature given by 
equation (fT7|) . As can be clearly seen, the steady state distribution is closely approximated 
by the Gaussian distribution fixed using the analytically-derived effective temperature and 
steady state average position coordinate. 
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FIG. 9: Steady state probability distribution P(x) for e = 0.3 and k = 0.1 23]. The horizontal 
coordinates are in units of the shift xq and with origin at the steady state average position xq{Pn+i)- 
The numerical solution is given by the solid line and the Gaussian fit is the dashed line. 

C. Mechanical resonator dynamics in the weak coupling regime 

In the previous section, we found that for weak coupling (n <C 1) the SET appears to 
the oscillator in the steady state as a thermal bath with effective temperature ^b^set = 
e Kis(-P)v(-P)A r +i- Suppose now that the resonator is not in a steady state, e.g., it is given 
some initial displacement amplitude and released, undergoing subsequent non-steady state 
motion. Does the SET still appear to the oscillator as a thermal bath? The following 
analysis addresses this question and as we shall learn, the SET indeed behaves as a bath, 
provided there is a wide separation in the oscillator and SET dynamics timescales (e < 1) 
in addition to weak coupling. 

We seek some approximate way to solve the dimensionless master equations (J7|) and (JHJ) 
which takes advantage of the weak coupling condition «: < 1. The solution should further- 
more be a good approximation for times much longer than the oscillator period 2ir/u in order 
to establish that the SET damps the oscillator. The obvious procedure is to solve for the 
evolving probability distribution perturbatively using k as the small expansion parameter. 
Let us first put the master equations in the following concise, matrix form: 



V = (H + V)V, 



(19) 
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where V 



P N (x,v,t) 
P N+ x(x,v,t) 




+ i ~Er- k(x) E L - n(x) j (2Q) 

e2x i; ~ v £ I V e r + K ( x ) ~ e l + «<*) 



Vi = kx | 1 1 | (21) 



and 




M, = ^ 1 o'w-i 1 ' (22) 

where we have redefined the position coordinate such that its origin coincides with the 
steady-state value (x) ~ Er [see equation Equation (JTHJ) resembles the time- dependent 

Schrodinger equation, but without the imaginary i since it describes a classical and not a 
quantum system. The 'Hamiltonian operator' T~C gives the free, decoupled evolution of the 
independent oscillator and SET systems, while the operator V = V\ + V2 describes the 
interaction between the two systems with V\ giving the dependence of the SET tunneling 
rates on the oscillator position and V2 giving the SET island number dependence of the 
electrostatic force acting on the oscillator. 

Given the resemblance of equation (fTT)|) to the Schrodinger equation, we can consider ap- 
plying approximation schemes that have been developed in quantum mechanics. We shall use 
a method developed for open quantum systems which is sometimes called the 'self-consistent 
Born approximation' (SCBA). A system is open when it is coupled to another system with 
an infinite number of degrees of freedom. The latter, infinite system is commonly termed the 
'environment' or 'reservoir' and the approximation method seeks to derive simpler, effective 
equations of motion for the finite system alone where the environment degrees of freedom 
have been integrated out. In our case, the finite system of interest is the mechanical res- 
onator modeled as a harmonic oscillator, while the environment comprises the tunneling 

n 

electrons in the SET. Applying the SCBA as described in section 3.1 of Ref. [44], we ob- 
tain the following approximate effective equation of motion for the probability distribution 
Piio(z,v,i) of the oscillator: 

P{x, v, £)ho = HroPho(x, v, t) 

21 



+e w Tr SE T [V(t)P SET (0)] e~ n ^ l Pno{x, v, t) 

dt' e WHO 'Tr SET [V(*)Pset(0)] Tr SET [V(O^set(0)] e'^ 1 P no {x,v,t) 



o 

+ / eft' e WHO *Tr S ET [V(t)V(t')V SET (0)} e'^Paoix, v, t), (23) 
Jo 

where Tiuo — e2a; ^; — is the Hamiltonian operator for the free harmonic oscillator and 
V(t) = e~ n<it V 'e +71 - ' is in the interaction picture. The initial, t — probability distribution is 

( Pn(0) 

assumed to be a product state: V(0) = Puo(x, v, O)'Pset(O), where "Pset(O) = 

\Pn + i(0) 

The influence of the SET on the oscillator is contained in the 'trace' terms, such as 
Tr SET [V(*)7W(0)] = ELi Ej=i V a/3 (t)P S ET/3(0). All such terms in © and their integrals 
with respect to t' can evaluated analytically. The resulting explicit expression for (J23j) is 
rather complicated, involving t-independent terms, t-dependent decaying terms of the form 
e _i , as well as decaying oscillatory terms of the form e _i cos(et) and e _t sin(et). The terms 
also depend on the initial state Pset(O) of the SET. However, if e 1, then we can neglect 
all decaying terms involving e~* since, on timescales of order the mechanical period and 
longer, t > e -1 , such terms have a negligible effect on the oscillator dynamics. Furthermore, 
the dependence on the initial state of the SET drops out. Neglecting the time- dependent 
decaying terms is usually called the 'Markov' approximation, while the combined steps of the 
SCBA and Markov approximation are often simply called the 'Born-Markov' approximation. 
As a result, we obtain the following much simpler effective equation of motion for the 
oscillator: 



OP 
~dt 



e2% TT - V TT + Ke V {v~x) + e A (P) N (P) N+1 — — + — 
ov ox ov ov \ ov ox 



P, (24) 



where (P) n and (P) n + i (=1— (P) n ~ Er) are the steady state SET island electron number 
probabilities [see equation (JT3J)] and we have dropped the 'HO' subscript on the oscillator 
probability distribution P(x,v,t) for notational convenience. Expressing (|24j) in terms of 
dimensionful coordinates, we obtain 



OP 
~dt 



U 2 R X- V— + 7SET7^-^ + 



3d d 7set^b7set d 2 



(25) 



dv dx dv m dv 2 

where recall the renormalized oscillator frequency is lor = y/l — k u, the damping rate is 

7set = Keu (26) 

and the effective SET temperature Tset is defined in (|T7|) . 
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Equation ([25)1 . which goes by the name 'Klein-Kramers' or 'Fokker-Planck' equation |45j | . 
describes the brownian motion of a harmonic oscillator interacting with a thermal bath; 
the oscillator experiences a damping force due to the SET [the third term on the right- 
hand-side of equation ()25))] with quality factor Qset = ^/7set = l/( Ke ) (since k C 1, the 
renormalized frequency can be replaced by the unrenormalized frequency in the definition 
of the quality factor) and an accompanying Gaussian distributed thermal fluctuating force 
[the fourth term on the right-hand-side of (j2*3jl . called the 'diffusion' term]. The d 2 P/dvdx 



term in equation (}2"2|) (called the 'anomalous diffusion' term in Ref. 44j) has been neglected 
because it is of order e smaller than the diffusion term when time is expressed in units of 
the oscillator period; it should have only a small effect on timescales of order the mechanical 
period or longer. The first moment of the fluctuating force which gives rise to the diffusion 
term vanishes, (F(t)) = 0, while its second moment is 

(F(t)F(t>)) = 2lS * TkBTswr 5(t-t>). (27) 
m 

The delta function in equation ()27|) signifies that the kicks inflicted on the oscillator by the 
SET are uncorrelated, no matter how small is the distinct time interval separating the kicks. 
This amounts to taking the limit T tunn ei — > 0, which is justified provided Ttunnei *C 2-k/u 
(equivalently e < 1) and we are interested only in the oscillator dynamics on timescales 
t > 2ir/u as discussed above. The condition e < 1 appears necessary for the oscillator to 
perceive the SET as a thermal bath; if the condition did not hold, then it would be possible 
to infer from the oscillator dynamics that it was coupled to a tunneling electron device and 
not to a many degree-of-freedom system in a state of thermal equilibrium. 

The second and third terms on the right-hand-side of equation (}23j) involving traces of 
single interaction potential operators V(t) vanish by virtue of the Markov approximation 
and also the fact that the position coordinate was originally redefined such that the origin 
coincides with the steady-state value. The damping and diffusion terms in equation ()24|) 
arise from the fourth term on the right-hand-side of equation (}2*3*|) involving the trace of 
the interaction potential operator product V(t)V(t'). More specifically, only products of the 
form V2(t)Vi(t') and V^i)"^^') [see definitions (|2"T|) and (}22|) ] give non-zero contributions, 
with the damping and frequency renormalization terms arising from the first product and the 
diffusion and anomalous diffusion terms arising from the second product. Thus, the effective 
thermal bath description for the SET crucially requires one to take into account not only 
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the influence of the SET on the oscillator (V2), but also the influence of the oscillator on the 
SET as well (Vi). 

In our analysis of the coupled SET-oscillator dynamics, we have neglected the coupling 
between the oscillator and its external thermal bath. Incorporating the external bath into 
the equations of motion is straightforward: we just add the term 

d(vP) lcxt k B T ext d 2 P 
7cxt^ 1 -^-r (28) 

to the above Fokker-Planck equation (J25j) . where we assume again for simplicity that the 
influence of the external bath on the oscillator is uncorrelated and memoryless. The damping 
and fluctuating force terms due to the SET and external bath can be combined respectively 
to yield a single damping term with effective quality factor = Qg^ T + Q~ x \ and a single 
fluctuating force term with effective temperature 

T eff = Q eff ( ^1 + . (29) 

VVSET Vext/ 

Summarizing so far, we have learned that provided the mechanical oscillator and SET 
are weakly coupled (k <C 1) and, furthermore, provided the SET dynamics occurs on much 
shorter timescales than the oscillator dynamics (e < 1), then the oscillator behaves effec- 
tively as if it is in contact with a thermal bath, with figure U\ being replaced by the much 
simpler figure ITU1 

For typical SET device parameters and micron-scale resonators with frequencies of the 
order of a few megahertz, separated by a resonator-island electrode gap of about 0.1 fim, one 
obtains effective SET temperatures of around IK and effective SET quality factors ranging 
from around 10 4 down to 10 2 as the gate voltage increases from 1 V up to 10 V 2^. Ex- 
ternal temperatures and quality factors for micron-scale resonators are of the order 100 mK 
and 10 4 -10 5 , respectively, so that it should be possible to measure the effects of the SET 
temperature and quality factor on the mechanical resonator. Experiments are underway 
in Keith Schwab's group to probe the coupled dynamics of the SET and mechanical res- 
onator 4f|. They have seen evidence in several devices that the SET behaves as an effective 
bath, heating the resonator and causing damping which increases with gate voltage. This 
is the first time that the back action of electronic shot noise on a nanomechanical resonator 
has been observed. Note, however, that they work with superconducting SETs (SSETs), 
where tunneling processes involve Cooper pairs as well, so that the above formulae for T SET 
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^eff ^SET ^ext 



T eff =Q eff (T SET /Q SET+ T x /Q ext ) 



Effective heat bath 



FIG. 10: Scheme of the SET-oscillator system under the conditions of weak coupling and 
wide separation of dynamic timescales between the SET and oscillator. The oscillator under- 
goes thermal brownian motion, behaving as if in contact with a thermal bath at temperature 
T e s = QcS (Tset/Qset + T ext /Q e xt), where Q~ s = Qset + QeA- 



and Qset do not app 
coupled dynamics (4' 



It is important then to also analyze the SSET-mechanical resonator 



IV. EFFECTIVE EQUILIBRIUM DYNAMICS OF OTHER NEMS 

In this section, we consider the effective equilibrium dynamics of some of the other NEM 
devices described in section [H] One of the first theoretical investigations to establish that 
a far-from-equilibrium electronic device can behave like an effective thermal bath was con- 
ducted by Dima Mozyrsky and Ivar Martin at Los Alamos 22]. They considered a model 
device comprising a single, electrical tunnel junction with tunneling rates depending on the 
position coordinate of a nearby mechanical oscillator, as illustrated in figure El Solving the 
quantum Schrodinger equation for the coupled, tunneling electrons-oscillator system, they 
found under the conditions of weak coupling and Markovian approximation that the oscil- 
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V 



FIG. 11: Oscillator coupled to an electrical tunnel junction [22j|. As a result of an applied voltage 
V, electrons will tunnel from the right (R) to left (L) electrode, with tunnel rates depending on 
the position a; of a nearby mechanical oscillator. The electromechanical coupling can be achieved 
by putting a net non-zero charge (or voltage) on a metallized mechanical oscillator. 

lator behaves effectively as if it is in contact with a thermal bath, with the tunnel junction 
characterized by a damping rate and a temperature given by 



B-i TJ 



eV/2, 



(30) 



where V is the voltage applied across the tunnel electrodes. 

Thus, despite the differences between the tunnel junction and SET (two tunnel junctions 
in series with gated central island), both can behave effectively as a heat bath. In fact, 
the similarity in their statistical properties goes further. It turns out that their respective 
temperature expressions (JT7jl and (}3Tjj) can be commonly equated to the ensemble-averaged 



energy gained by an electron due to tunneling across a junction 



23j. For the tunnel junction, 



the averaged energy gained is eV/2, assuming zero electron temperature and constant density 



of states in the electrodes, as well as energy-independent tunneling matrix [22, For 
the SET, with the same asumptions the average energy gained by a tunneling electron is 
(E R /2)(P) N + (E L /2)(P) N+1 « eVd s (P) n{P) N+i for k <C 1, where we have used equations 
©, (HOD and (H3J). 

While the mechanically compliant tunnel electrode illustrated in figure ^ closely resem- 
bles the just-considered model tunnel junction (figure 1X1)1 . there is a difference in that the 
tunneling electrons will give a momentum recoil to the mechanically compliant electrode, 
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which is of course absent for the fixed electrode. To our knowledge, there has yet to be 
a comprehensive analysis of the coupled dynamics of tunneling electrons and mechanically 
compliant electrode, including the effects of momentum recoil. Nevertheless, one might 
expect that under the appropriate conditions of weak coupling and wide separation of dy- 
namical timescales between the mechanical and electronic degrees of freedom, the effective 
thermal bath description again arises. 

As mentioned in section |HJ the electronic and mechanical dynamics are strongly coupled 
in the classical charge shuttle (figure in order for a current to flow, the nanomechanical 
island electrode must shuttle charge between the drain and source electrodes, so that the 
dynamical timescales of the mechanical and electronic degrees of freedom are comparable. 
Because it is not possible to have a wide separation of timescales, it is unlikely that there is 
a regime in which the electron tunneling dynamics is effectively that of a thermal bath, as 
perceived by the mechanical shuttle. The only possibility of separating the timescales is to 
make the drain-source electrode gap smaller than the electron tunneling length (so that a 
tunneling current can flow without the nanomechanical island electrode having to move from 
drain to source electrode), likely an unrealisable regime for nanoscale mechanical resonators 
that must fit in the gap. 

Thus, provided the electron dynamics is stochastic, the electromechanical coupling 
strength is weak and the mechanical dynamics is much slower than the electron dynam- 
ics, the above examples suggest that the electronic device can be effectively replaced by a 
thermal bath for the reduced dynamics of the coupled mechanical resonator, so that fig- 
ure [3 can be replaced by figure El (but now refering to a generic NEMS). How might one 
go about providing a proof of this conjecture for a general class of electromechanical sys- 
tems? A possible direction is suggested by the Born-Markov approximation method applied 
to the SET- resonator system in the previous section IIII CI In fact, the derivation of the 
Fokker-Planck equation (that ensures the effective thermal bath description) does not de- 
pend on the detailed form of the matrices appearing in the definitions (J2T|) and (j2Zj) of 
the interaction operators Vi and V2, respectively; the derivation of the Fokker-Planck equa- 
tion should carry through for additional and different types of tunneling process, e.g., for 
a superconducting SET 4s|. Of course, the detailed expressions for the renormalized fre- 
quency, damping constant and temperature will differ from one device to another. The 
important necessary steps that lead to the Fokker-Planck equation are the second order ex- 
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pansion in the interaction operator V (Born approximation-requires weak coupling) and the 
neglect of time-dependent decaying terms in the electron dynamics (Markov approximation- 
requires fast electron dynamics as compared with mechanical oscillator period). A suitable 
starting point for a general proof would be to write down a Boltzmann/master equation 
involving a probability density function Pi^\(x,v,t) where (x,v) are the oscillator's clas- 
sical position/velocity coordinates and {k} denotes some appropriate choice of electronic 
coordinates. The general interaction operator V should involve both position and velocity 
(to account for momentum recoil) dependent terms. The proof would then involve apply- 
ing the Born-Markov approximation to this master equation, recovering the Fokker-Planck 
equation. An important proviso, however, is that there is no guarantee that the electronic 
device can necessarily damp the oscillator. The tunneling dynamics may be such that it is 
more likely for an electron to give rather than to take energy from an oscillator, resulting 
in unstable coupled dynamics. This instability manifests itself through 'negative damping' 
and 'negative temperature' coefficients in the Fokker-Planck equation. Such an instability 
occurs in the case of the superconducting SET for certain ranges of gate and drain-source 
voltage 42, 3- 

A proof of the effective thermal bath description along; the above lines has recently been 
demonstrated by Aashish Clerk at McGill University 40(. He gives a quantum treatment, 
with the electronic device modeled as a generic linear response detector of the quantum 
oscillator's position (essentially the quantum version of figure EJ), to which it is weakly cou- 
pled. The proof is in part a generalization of Mozyrsky and Martin's analysis of the tunnel 
junction |22j. However, as described above, when the coupled dynamics under consideration 
is classical it should be possible to give a similar proof within the framework of classical 
master equations. Indeed, in Mozyrsky and Martin's tunnel junction analysis, the condi- 
tion eV 3> ftw is assumed, which from equation ()3U|) amounts to taking the classical limit 
&bTtj 3> hu for the oscillator. Certainly, a classical derivation of the effective thermal bath 
description for NEMS avoids having to simultaneously deal with the complicating (but fas- 
cinating) issue of how classical dynamics arises by approximation from quantum dynamics. 



28 



V. CONCLUSIONS 



We have presented a brief overview of the field of NEMS research, with an emphasis on 
the classical effective dynamics of a nanomechanical resonator due to its interaction with 
a tunneling electron system. Under the conditions of weak coupling and the characteristic 
mechanical dynamics occuring on much longer timescales than the tunneling electron dy- 
namics, the electron system can behave effectively as a heat bath, causing the resonator to 
undergo thermal brownian motion. 

As they have been stated, the conditions for the resonator to perceive the electronic 
system as a thermal bath are rather general. It is natural therefore to wonder whether 
there are other (semi) classical examples involving (not necessarily nanoscale) mechanical 
oscillators coupled to (not necessarily electronic) far from equilibrium systems with many 
degrees of freedom, such that the oscillator undergoes effective thermal brownian motion. 
Many such examples have in fact been considered. One example is laser Doppler cooling of 
a harmonically trapped ion The laser is tuned below the atomic resonance, so that 

the Doppler shift makes it more likely for the atom to absorb a photon with momentum 
oppositely directed to that of the atom than in the same direction, resulting in damping of 
the atom's motion. Despite the fact that the laser radiation field is obviously not a thermal 
equilibrium state, nevertheless under certain weak coupling and wide timescale separation 
conditions, the atom's dynamics is described by a Fokker-Planck equation j^l|, with effective 
temperature much lower than the ambient temperature: the laser cools the trapped ion. 
Another, recently considered macroscopic oscillator example j^] involves placing a sphere 
(ping-pong ball) in an upward gas-flow. It was found that the sphere behaves effectively as a 
harmonic oscillator undergoing thermal brownian motion, despite the far-from-equilibrium, 
turbulent state of the gas molecules. A derivation of this behavior might start with a 
Boltzmann equation for a single massive particle undergoing collisions with a dilute gas 
of small mass particles, and then identifying the relevant small dimensionless coupling and 
timescale parameters for a subsequent series expansion of the Boltzmann equation. However, 
the recovery of a Fokker-Planck equation is not straightforward, owing to the relevance of 
the non-trivial turbulent gas dynamics for the effective dynamics of the oscillator. Yet 
another, related macroscopic example involves immer sing a torsion oscillator in a far-from- 
equilibrium, vibration-fluidized granular medium Again, the granular medium 
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was found to behave effectively as a thermal bath, with the oscillator undergoing thermal 
brownian motion. 

Returning to the subject of NEMS, there are many aspects of their dynamical behaviour 
still to be understood. In the strong coupling regime (e.g., k > 1 in the case of the SET- 
oscillator system), little is known in general about the effective dynamics of the oscillator 
and also how the coupled dynamics manifests itself in the effective behavior of the electronic 
system (e.g., current, current noise, etc.). Furthermore, little is known in general about 
the quantum dynamics of NEMS and how (semi) classical dynamics arises as a result of 
dephasing within the electronic system, as well as due to the external environment of the 
nanomechanical resonator. 
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